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Abstract 

Let A be the finite difference Laplacian associated to the lattice Z d . For dimension d > 3, 
a > and L a sufficiently large positive dyadic integer, we prove that the integral kernel of 
the resolvent G a := (a — A) -1 can be decomposed as an infinite sum of positive semi-definite 
functions V n of finite range, V n {x—y) — for \x— y\ > 0(L) n . Equivalcntly, the Gaussian process 
on the lattice with covariance G a admits a decomposition into independent Gaussian processes 
with finite range covariances. For a = 0, V n has a limiting scaling form L^ n ^ d ^ 2 ^T c ^ (^rr) 
as n — > oo. As a corollary, such decompositions also exist for fractional powers (— A)~ Q / 2 , 
< a < 2. 

The results of this paper give an alternative to the block spin renormalization group on the 
lattice. 

keywords: gaussian processes, finite range decomposition, lattice, renormalization group, Levy 
processes. 

1 Introduction 

A smooth Gaussian process ((x) on R d with the property that the expectation E£(x)£(y) = 
when \x — y\ > L will be said to have finite range L. What is the class of Gaussian processes (f> 
that can be expressed as a sum <f> = Y2j Cj °f independent finite range processes with ranges ~ V 
for some L? Let us call such processes finite range decomposable. 

We can reformulate this in terms of the covariance: a Gaussian process cf) is finite range decom- 
posable if the covariance C{x,y) := Ecfi(x)(fi(y) can be written as a sum C = ^2jVj where each 
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Vj(x,y) is positive semi-definite and has finite range ~ U . In this form the question has already 
received a partial answer in the study of ground states for many-body Hamiltonians. In particular, 
in |HS02| Hainzl and Seiringer discuss this background and consider the decomposition 

POD 

V(x)= dr g(r)x r /2* Xr/2(x) (1-1) 
J 

of a radial function V{x) as a weighted integral of tent functions x r /2 * Xr/2{ x )i where Xr/2 ls the 
indicator function of the ball of radius r/2. An explicit formula for g in terms of V is derived. For 
example, in three dimensions, 

g(r) = _ V"(r)A0' 

7T 

so necessary and sufficient conditions for g > in terms of V are readily formulated. In particular 
Coulomb and Yukawa potentials in three dimensions have decompositions with nonnegative g. 

This is relevant to our question because the tent function is positive semi-definite and therefore, 
when g(r) > and I is an interval [a, b), 

Vi(x) := j^dr g{r)x r /2 * Xr/2 0) 

is also positive semi-definite. By breaking up the range of the r integration in Ijl.lJ) into a disjoint 
union of intervals, Ij := [LP, LP + ),j £ Z, we have V = Yl Vl~ an d there is a corresponding finite 
range decomposition <j) = J2j Cj when <f> is the Gaussian process with covariance V{x — y) with 
g(r) > and Cj has covariance Vj.. 

These decompositions are not the final answer to our question, because we are also interested 
in kernels defined on the lattice Z d and furthermore one may get a wider class by not insisting 
on decompositions based on tent functions. For lattices or the continuum we have preliminary 
results that suggest that resolvents of quite general elliptic operators and fractional inverse powers 
of elliptic operators are candidates for such decompositions. 

Our interest in this question is rooted in the Renormalization Group (RG). In quantum field 
theory and other contexts the RG is a method to calculate the expectation EZ of a functional 
Z = Z(<p) of a Gaussian field </>. One decomposes <f> = Ylj>i Cj as a sum °f independent Gaussian 
fields Cj and integrates out each (j one at a time. Let Ej be the expectation that integrates out Cj- 
Then the RG is the sequence of maps Zj t— > -Zj+i := Ej + \Zj. EZ is obtained from EZ = limZj, 
starting with Zq := Z. [These ideas are explained further in Section E]. The point is to choose the 
decomposition to have special properties so that each expectation Ej is more amenable to analysis 
than the whole expectation E and furthermore so that the map Zj t— > Zj + \ can be analysed within 
the context of dynamical systems. In particular the RG is very informative when the limiting map is 
autonomous, up to a scaling. This is possible when the covariances Vj of Cj are becoming self-similar. 
This means that there should exist a dimension [(f)] such that Vj(x,y) =: L^ 2 ^Tj{L^^x,L^^y) 
defines a scaled covariance Tj which tends to a limit as j — ► oo. 

When the covariances Vj are finite range, the map Zj t— > Zj+i can be studied by using the 
independence of Cj( x ) an d Cj(u) f° r \ x ~ u\ > L 3 . We amplify on this remark at the end of this 
introduction. In some ways, the use of these finite range covariances gives the simplest framework 
that goes beyond the hierarchical models. This program is also close to the technique that was 
invented by Frohlich and Spencer in |FS81j to study the Kosterlitz-Thouless transition, but it is in 
principle more precise and more robust. 
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Let A be the finite difference Laplacian associated to the lattice Z d . We will consider decom- 
positions for the kernel G a (x — y) of the resolvent (a — A) -1 and we will also consider the Green's 
function of a stable Levy process which is the kernel C(x — y) of (— A)~ Q / 2 where < a < 2. The 
mass parameter a £ [0, oo). We state our results for dimensions d > 3 because lower dimensions 
require extra discussions for the case a = 0, but the basic construction is valid in lower dimensions 
as well. L is a parameter of the form 2 P . p can be any sufficiently large integer. There are many 
technical results in this paper so we have summarised the main points in the following theorem 
which is a combination of results from the theorems in the rest of the paper. 

Theorem 1.1. For n = 0,1,2 ... , and all a > 0, there are positive semi-definite functions r^(x) 
defined for x G (X~ n Z) such that 

1. G a {x -y) = J2 L- 2n[<t>] rf " a with [4>] = {d- 2)/2 

2. T a n (x) = for \x\ > 6L 

3. \f a n {p)\ < c KL {l + \p\)~ k forpe [-L n n,L n 7r] d 
4- i"c*(p) := lhrin-^oo ^n(p) exists pointwise in p 

5. Fix a positive integer I and let e = L~ l . Then = linin^oo r° exists in L°°((eZ) d ). Further- 
more the L°° limit of any multiple lattice derivative ofT^ also exists and is the corresponding 
continuum derivative ofT^ 

6. Analogous statements hold for C{x — y), but with [</>] = (d — a)/2. 

The multiscale expansion in (1) above is obtained in (j3.31|) of section 3, and the finite range 
property (2) is given in Lemmas 13.11 13.21 The bound in (3) follows from Theorem 15.51 of section 
5. (4) above is obtained in section 6 in the course of proving Theorem 16.11 see ([6.24)1 . (5) above 
is part of Theorem 16.11 Turning to (6) above, the finite range decomposition of C is obtained in 
section 4 ( see (|4.2|) . 1)4. 3 jl et seq.). A uniform bound on Levy fluctuation covariances is supplied 
in section 5 ( see Corollary 15.61 ). The statement analogous to (5) above is part of Corollary 16.21 of 
section 6. 

Let us call the Gaussian fields in the decomposition fluctuation fields. Other (wavelet) decom- 
positions were developed in the context of the Block Spin Renormalization Group of Kadanoff 
and Wilson by Gawedzki and Kupiainen, |GK80l IGK83j as well as, in related work, by Balaban, 
|Bal82al IBal82bj . Although the Gawedzki-Kupiainen fluctuation fields are not finite range, they 
have their own advantages: notably they are independent lattice fields, determined by random vari- 
ables defined on increasingly coarse lattices of spacing L n . Our decomposition achieves this only 
in the weaker sense that, with high probability, Cn( x ) ~ Cn{y) f° r \ x ~ y\ L n , but £ n retains low 
probability variations on all smaller scales. This is a price paid for retaining translation invariance 
on small scales. 

In section 3 the aforementioned finite range decompositions are obtained. The rescaled fluctu- 
ation covariances live on finer and finer lattices, but all have the same finite range. In section 4 
probabilistic aspects of our construction are discussed and the finite range decomposition of the 
Levy Greens function is obtained. We are interested in the Levy Greens function (— A)~ Q//2 be- 
cause varying the parameter a affects the scaling of the associated field and gives insight into the 
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dynamical system Zj \— > 2>j+i. Moreover it is of intrinsic interest in various problems in probability 
theory. For example we may wish to study critical properties of self avoiding Levy walks by renor- 
malization group methods. We also discuss in this section renormalization group transformations 
based on the above finite range decompositions. 

Section 5 is devoted to bounds. Here our main result is Theorem 15.51 This theorem states 
that every member of the sequence of rescaled fluctuation covariances is uniformly bounded in 
lattice Sobolev norms of arbitrarily high degree. The bound is independent of the lattice spacing. 
Finally in Section 6 we prove that the sequence converges in Sobolev norms in a precise sense to 
its continuum limit which is appropriately identified. The continuum limit is smooth. This is the 
content of our main Theorem 16.11 

We conclude with a brief indication of the role of the finite range property in the analysis of 
Zj i — ^ -Zj+i- Consider Zq \— > Z\ when 

Z = Z (A,fo) :=n^ W 

where A C Z d is a large box shaped subset of lattice points which is a disjoint union of some 
standard cube shaped subsets of lattice points A C Z d of side greater than the range of the first 
field Ci m a finite range decomposition </>o = Ylj>i Cr Let (\>\ := YLj>2 0> X C Z d and 

Zi(X) = Yl e - A ^M 
xex 

Note that Z\{X) is independent of Q\. Then 

Zx(A) = E x J] Z (A) = Ex J] (z (A) - ^i(A) + Z X {A)\ 

AcA AcA ^ ' 

Write 5Zq (A) := Zq(A) — Zi(A). The product expands into a sum over X of terms 

Zx(A) H 5Z (A) 

A<£X AcX 

In other words, X labels the factors where 5Zq is selected. We can partition X into disjoint 
connected components Xx, ■ ■ -Xm, where Xi is connected when the cubes A in Xj are such that 
one can pass between any pair of cubes by a path whose steps are nearest neighbour cubes in Xi. 
Let 

K (Xi) := 11 SZ (A) 
AcX, 

Then 

2i(A) = J]^ 2 1 (Xo)E 1 l[K (X i ) 

' Xi,...,X M i 

where the connected sets Xj are disjoint and Xq := A \ Uj>iXj. Notice that Z\ can be and has 
been moved outside the expectation since it is independent of (j. Now comes the key point: The 
sets Xj are connected unions of nearest neighbour cubes of side length greater than the range of 
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Ci- Since they are disjoint they are separated by a distance greater than the range of £i- Therefore, 
by the finite range property, 

2i(A) = J]^ Z 1 (X Q )l[E 1 K (X i ) (1.2) 

' Xi,...,X M i 

The perturbation on a large volume A has been reduced to local calculations E\KQ{Xi) and the 
standard but heavy machinery of cluster expansions is being replaced by independence and geom- 
etry. 

Unlike Zq(A) the image functional Zi(A) no longer factors into contributions from boxes, which 
is the great simplification of hierarchical models, but there is still a large part Xq of A where this 
property is retained. For the next RG map one proves that the more general form (|1,2|) is stable 
in the sense that Zi can also be written in the same form but with a different K\ and with larger 
cubes LA. The program is then to prove that EjKj(X) gives very little weight to connected sets 
X which are unions of many boxes LP A. This can be facilitated by making a better choice of Z\, 
since the derivation is valid for other choices of Z\. In particular, one can replace A in Z\ by some 
other value Ai chosen to minimise E\Kq. This idea leads to a flow of the coupling constant A — > Ax- 

The first use of finite range covariances was in |MS00j . |BMSj is another appearance of finite 
range covariances. 



2 Preliminaries 

Throughout we will assume that d > 3. Let L be a large integer power of 2. Define e n = L~ n . We 
will be working on a sequence of lattices (e n Z) d c R d , {e n Z) d C (e n+1 Z) d , with n = 0, 1, 2, ... and 
eventually passing to H d . (e n Z) d is equipped with the discrete topology. The measurable sets are 
subsets of points and the measure dz on (e n Z) d is defined by 

/ dz f(z) = e d n /(*) (2- 1 ) 

We endow H d with the distance function 

\x-y\ = maxi<j< d \xj - yj\ (2.2) 

Let 

TJ(-P\ — ( 

2 ' 2 



U(R) = (-^,^) d CK d (2.3) 



be an open cube of edge length R. Define 

U En (R) = U(R) n(e n Z) d (2.4) 

and its boundary 

du e n = {y U En : \x-y\= E n , some x € U £n ] (2.5) 

The distance function l.l is that induced from H d . We denote by U F = U F U dU F the closure of 
U £n . The lattice Laplacian A En is defined by the quadratic form 

(/, -A £n /) i2((£ftZ)d) =E d n J2 £ n 2 \m - f(y)\ 2 (2-6) 

<x,y> 
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where the sum runs over the nearest neighbour points in (e n 7i) d . Let S = {e\, ...,id} be the set 
of standard unit basis vectors in Z d . For any such lattice unit vector e E S define forward and 
backward lattice derivatives V±e of a function in the direction e by 

(V± s /)(x) =e-\f(x±e n e)- f{x)). (2.7) 

The backward derivative is defined so as to be the adjoint of the forward derivative. Then the 
definition ()2.6[) for the lattice Laplacian can be written as 

(/,-A en /) i2((£nZ)d) = 4 J2 K V eiOMI 2 ( 2 -8) 

The corresponding resolvent G £n with a > is 

G a £n {x-y) = (-^e n + a)-\x-y) (2.9) 

P 6 (2.10) 



[-7r/e n ,n/s„] d ( 27r ) d a-A £n ( P ) 

where 

d 

A £n (p) = 2e~ 2 £ (cos(e„^) - l) (2.11) 

At=l 



3 Multiscale Decomposition of the Resolvent 

We say that a function f(x,y) has finite range R if 

f(x, y) = for \x — y\ > R. 

Consider the resolvent in Z d 

G a (x-y) = (-A + a)- 1 (x-y) (3.1) 

with a > 0. We will first develop a multiscale decomposition for the resolvent G a , a > into 
smooth finite range positive semi-definite functions. 

As in l(2~3|) and U(L) C R ' is an open cube of edge length L in K d and U E (L) = 

U(L) n (eZ) d the induced cube in (eZ) d with dU e {L) its boundary. Assume that the cube is 
centered at the origin. We will suppress the argument L when there is no risk of confusion. 

On the lattice there is no need to distinguish functions from measures. Nevertheless we use 
measures in cases where the associated continuum object is a measure. A case in point is the lattice 
Poisson kernel Vjj s (x,du), which by definition is the measure supported on dU e such that 

h(x) = (x, /) := J V% e (x, du)f(u) (3.2) 

is the unique solution to the boundary value problem 

(-A e + a)h(x) = : x £ U e (3.3) 
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h(x) = f{x) : x G dU £ 



(3.4) 



where / : dU £ — ► R. Existence and uniqueness are easily proved since h solves a finite dimensional 
set of linear equations. Note that because a > a solution /i(x) satisfies the weak maximum 
principle. In Section 0] we will see an explicit construction of which shows that Vjj is a defective 
probability measure. Defective means that the mass is at most one. (For a = and / = 1, h = 1 
which implies Vjj {x, 1) is a probability measure). We will say h is a-harmonic in X £ if /i solves 
(— A e + a)h{x) = in X n (eZ) d . 

Let g(x) be a rotationally invariant non-negative C°°(R d ) function of compact support such 
that 

9(x) = : \x\ > - (3.5) 

with the normalization 

/ 5 (x) = 1 (3.6) 
We restrict g(x) to the lattice (eZi) d and choose the normalization constant c £ such that 

/ dx c £ g(x) = 1 (3.7) 

J{eZ) d 

Since g(x) is a continuous function of compact support Riemann sums converge. Hence c e is a 
continuous function of e on the compact set < e < 1 and thus uniformly bounded. Moreover 
c e — > 1 as e — > 0. 

Now comes the main idea. The point of the function g is to avoid needing detailed knowledge 
of the Poisson kernel for the lattice. We are about to use the Poisson kernel to define an averaging 
operator that leaves a-harmonic functions unchanged. This property leads to the finite range 
property in Lemma 13.11 It is relatively easy to prove that our averaging operator is smoothing 
(uniformly in lattice spacing) because when checking differentiability, derivatives either fall on g 
which is smooth by choice or on the x argument of Vjj (x, du) with x forced to be away from the 
boundary dU £ so that the easy part of standard elliptic techniques is sufficient to prove smoothness 
uniformly in the lattice spacing. 

Given a function / : (eZ) rf -tRwe define the averaging map : 

f - A a £ (L)f 

where 

{A a £ {L)f){x)= [ dzc £ g(z-x) I ' V^ {L) (x - z,du)f(u + z) (3.8) 
Note that this can also be written as 

(A a £ (L)f)(x) = [ dz c £ g(z - x)PZ, L) (x, f) (3.9) 

where U £ (L, z) is the translate of the cube U £ (L) so that its center is now z . In summing over the 
translates we have put in the smooth function g , and not the delta function, because we will need 
to take derivatives with respect to x and this is hard to do if x is the center of the cube. 
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Now this integration over all translates makes A £ (L) translation invariant. Translation invari- 
ance plays an essential role in the proof of positive semi-definiteness of the fluctuation covariance 
constructed below, (see Lemma f 3 . 1 1) . Proof: For b E (eZ) rf , let fb{x) = f(x — b). By change of 
variables z — » z — b in (|3,8j) 

(A a £ (L)f)(x-b) = (A a £ (L)f b )(x). (3.10) 

^From (j3.8j) we see that for every fixed x, (A £ (L)f)(x) defines a bounded, positive linear functional 
on Co((eZ) rf ), the space of functions of compact support on (eZ) d . We have 

\(A a e (L)f)(x)\ < H/IU 

so that the norm of this linear functional is < 1. This gives a family of defective probability 
measures A £ (L)(x, du) on (eZ) d : 



(A«(L)f)(x)= / A*(L)(x,du)f(u) (3.11) 

J(eZ) d 



The Fourier transform of this measure 



A a £ (p) = duA a E (L)(0,du)e- ip - u 

>(eZ) d 



satisfies 

\A a e (p)\ < f duA a e {L){Q,du) < 1 (3.12) 

J(eZ) d 

Define the fluctuation covariance 

T a e (x -y)= G a e (x -y)- (A a £ (L)G a £ A a £ (L)*)(x - y) (3.13) 

where by definition 

{A a £ {L)G a £ A a £ {LY)(x-y)= [ [ A a £ (L)(x,du)G a (u - v)A a £ (L)(y,dv) (3.14) 



The latter is the analogue of the block spin covariance in statistical mechanics, jGK86j . 

Lemma 3.1. T £ and A £ (L)G £ A £ (L)* are positive semi- definite. T £ has finite range, 

T a £ {x-y) = : \x - y\ > 3L (3.15) 

and the Fourier transform T £ (p) is continuous in p including at p = 0, uniformly in e. 

Proof. First we prove the finite range property. By the definition of the Poisson kernel, if / is 
a-harmonic in x + U e , then V^j (x, f) = f(x). Since c e g in the definition of A £ (L) was chosen to be 
a probability density with support in U e (L/4), for / a-harmonic in x + U e (5L/4), 



(A a £ (L)f)(x) = / dz c e g(z - x)V a Ve(L Ax, f) = dz c £ g(z - x)f(x) = f(x) 

J(eZ) d K ' ; J(sZ) d 

When \x — y\ > 3L, x + U £ {bL/A) and y + U £ (5L/4) are disjoint. Therefore G £ (u — v) is a-harmonic 
in each argument in the appropriate region and therefore 

(A a £ (L)G a £ A a £ (LY)(x-y) = G a £ (x-y) 
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which proves (|3.15j) . Now we prove positive definiteness. By translation invariance, we can take 
the fourier transform of T £ (x — y) to get 

t a £ (p) = (l-\A a £ (p)\ 2 )G a £ (p) (3.16) 

Now G £ (p) > and |A"(p)| < 1 where we have used (|3.12j) . Hence T £ (p) > 0. This proves positive 
definiteness of T £ . The positive definiteness of A £ G £ A £ * is obvious. 
Continuity of f«(p): By lUTl . 

G»(p) = 3 

a-A £n (p) 

is continuous if a > or p ^ 0. If a = then (|3.16j) shows that the p~ 2 singularity of T £ (p) * s 
cancelled by 

\l-A a £ (p)\ = \A a M-A a £ (p)\=o(p 2 ) 

which holds because p derivatives of A £ (p) are moments for A £ and this is a probability measure 
of compact support which therefore has moments of all orders. □ 

Choose e = e n -\ for n > 1 and write (|3.13|) in a rescaled form. It is easy to check that for 
u, v £ (en-iZ)^ we have 

G^_ 1 (n-,)=L-( d - 2 )G^(^) (3.17) 

^_w*.*«)=<Wz't) (3 - 18) 

Define the sequence of functions g n on R d by 

5n (z)=L" d <?(L"z) (3.19) 

where g is the function introduced earlier (see (|3.5|) . Q3.6|) ) and observe that because of the normal- 
ization H3.6|) . the function g n is also normalized : 

/ dz g n (z) = l (3.20) 

7Rd 



Moreover, from the support property of g, we have that 

4L 



9n(x) = : \x\ > — ^ (3.21) 



Let 

R m = L"^- 1 ) (3.22) 
As in Q3.11JI we have in (e n Z) d for n > m > the measure A" m (R m )(x, du) given by 

/ duA a £n:m (R m )(x,du)f(u)=[ dzc £n _ m g m {x-z)Vfj e{Rm ^{xJ) (3.23) 

Note that 

4Lo = < (3-24) 



9 



as defined earlier. Observe that 



dz c £n _ m g m (z) = dz c £n _ m g(z) = 1 (3.25) 

{e n Z) d J(e n - m Z) d 

from the definition of the constants c Sn in ((3.7|) , The definition (|3.19|) together with (|3.23f) and 
((3.18(1 imply the scaling relation 



A a en _ 1>m -i(Rm-i)(x,du) = A^ m (R m )(^) ( :j.2(i) 



Applying this to the righthand side of 1)3.13(1 we get for n > 1 

^e n -i (* - y) = O-i 0= " » ) + L_( ^ 2) OCiW < a Of) (^) (3-27) 
We can now iterate 1)3.27(1 starting with n = 1, n-times using the same principle. Define for n > 1 

n 

^ = II ^.n+l-i^^ ) (3-28) 
i=i 



For n = we set 

We also define on the (s n Z) d lattice 
and 



,48 = 1 

-pa _ jul ya a a* (o nq\ 



^ = ^G e a n ^r (3-30) 
Then we have the multiscale decomposition for the resolvent 

n-l 

G*(* - y) = L- jM if'"^) + L- n{d ~ 2) Sn" a (^) (3-31) 
3=0 

which is valid for a > 0. The special case o = gives the multiscale decomposition for the massless 
Green's function 

n-l 

G\x - y) = £ L-^~ 2 > + Ql (3.32) 

3=0 

Lemma 3.2. For all n > 

K(x-y) =0 : |s-y| > 6L 

Proof, is a multiple convolution, so the range of is the sum of the ranges of the convolved func- 
tions. ^From the definition ()3.23() of A^ n m (R), the support property ((3.21(1 of g m and (R)(x, u) 
vanishes if \u — z\ > R, we find that the range of A* n>m (R) is R+ jj^-t- From the definition (THHll 
of A% the range of A%_ is Y^=i L~( n ~i>(l + which is less than 3 for L large. By construction 
the range of T" is less than 3L by Lemma 13.11 The Lemma follows. □ 
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4 Probabilistic Aspects 



Multiscale decomposition for the Levy Green's function. 

Let xf^ , < a < 2, be the stable Levy process in Z d , (stable in the sense that its scaling limit 
is stable). Note that i^dxji |) < oo, provided a > 1. However in the following we will allow for 
the full admissable range < a < 2. The Levy Green's function C is given by 

C(x-y) = (-A)- a / 2 (x-y) = const da a - a/2 G a {x-y) (4.1) 



with < a < 2. In this range it is easy to verify that the integral representation converges. 

C has a finite range multiscale decomposition, which we obtain by inserting the multiscale 
decomposition 1)3. 31j) for G a into the integral representation 1)4. 1JI of the Levy Green's function C 
to get 

C( x -y) = Y J L ~ j(d ~ 2) const / da a" a/2 rp°(^) + 

+L -n(d-2) CQnst [°° da a -*/2 gLP-a t^pL} 

Jo L n 

After rescaling in a in each term we get 

71-1 _ _ 

C(x -y) = J2 L- 2M T A^r) + C n £j/) (4.2) 

3=0 



where 



oo 

roust I da a- a/2 T^ 





= const 




= const 


[4>] 


d — a 


2 



(4.3) 

Note that [(f)] as defined above is the canonical dimension of the the scalar Gaussian field <p dis- 
tributed with covariance C. We have [</>] > 0, since d > 3 and < a < 2. Tj and C n are well defined 
because of the bounds provided below (see section 5, Corollary I5.6|) . Moreover by Lemma I3~2l and 

L„(x - y) = : \x-y\>6L 

^From Lemma 13.21 and (|4.3|) . C n and T n are positive semi-definite and thus qualify as covariances 
of Gaussian measures denoted Ho n i^V n - The multiscale decomposition (|4.2|) now gives rise to 
renormalization group transformations. From Ij4.2|) we get for x, y S (e n Z) rf 

C n (x - y) = T n (x - y) + L" 2M C n+1 (^) (4.4) 
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and hence we have a sequence of RG transformations 

z n+l (4>) = / dn Tn (Q z n (C + <f> L -i) (4.5) 



where 

L - 1 (x)=L-M0(|) 

and 

J d Vc n ((t>) z n {4>) = j dii Cn+l {4>) z n+ i{4>) (4.6) 

with C = C given by (j4.1|) . 
Poisson Kernel. 

Let a4™^, i > be continuous time simple random walk with right continuous paths and state 
space (e n Z) d . The characteristic function is 

(n) 

and the generator of the Markov process x t is the lattice Laplacian A £n . Note that the semigroup 
e *A en j g a con traction on L°°((e n Z) d ). In the discrete topology the latter coincides with the space 
of bounded continuous functions. Hence e tAen is a Feller semigoup so that xf 1 ^ is strong Markov 
with respect to stopping times r. 

Now let x\ be the above process in (eZ) d . Let P x the probability measure for the process 
conditioned to start at x. As in ()2.3|) and (|2.4|) . U(R) is an open cube of radius R in R d , U £ {R) 
the induced cube in (eZ) rf , and the boundary dU e {R) is defined as in (|2.5jl . Let T[/ e be the first 
exit time from £7 £ . Then tjj £ is also the first hitting time of dU e from the interior. E x (tu s ) < oo 
because U e {R) is bounded, and hence tu £ < oo, P x a.s.. 

Let A be a subset of points in dU e . It is a standard result in probability that the measure 
Vjj e (x, dy) defined on dU e by 

V^(x,A) = E x (e- aT ^l XTUeeA ) (4.7) 

is the Poisson kernel we defined in (|3,2j) . 

Also, since P x {x Tu ^ S dU) = 1 the total mass is 

V?j e (x, 8U £ ) = E x (e- aT ^ ) < 1 (4.8) 

The same construction works in TL d with x\ replaced by standard Brownian motion xt, U £ {R) 
replaced by U(R) and A taken to be any Borel subset of dU(R). 

Let now R = R m = L^™ -1 ), with < m < n. In Section 6 we will need an estimate for 
1 — Vijr (x,dU en ) where it is understood U Sn = U £n {R m ). Observe that 1 — e~ aT < ar so 

0<l-V% en (x,dU En )<aE x (T Ue ). 

We can estimate the mean exit time as follows. Let / be a smooth function in R d bounded in a 
neighbourhood of U(R). Then by the strong Markov property, for any x G U £n (R m ), 

E x (( TUe dt(A £ J)(x t )) = E x (f(x TU )) - f{x) 
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We choose x = and f(x) = \x\ 2 in U(^R m ). Then by a simple computation in U(R m ) we obtain 
A £n f(x) = 2d. Moreover as shown in Section 6, (|6.14j) . dU Sn {R m ) C dU(R). Hence we get the 
bound 

E x (r Ue ) < ±E x (xl r ) < ^. 



Thus we have proved the following 

Lemma 4.1. Let R m = L~( m ~ i \0 < m <n Then we have the bound 

0<l-^(K m )(^^ £n ( J R m ))<a^ 

Now recall the definition of the measure m (i? m )(0, du) on (e n Z) d given in (|3,23|) . c £n _ m g m (z) 
is a probability density in (eZ) d . Hence from (Lemma 14. If) we get 

Corollary 4.2. 

0<1- / A« n>m (R m )(0,du)<a^ 



5 Bounds 



We will give uniform bounds on Fourier transforms and Sobolev norms of arbitrary high index. In 
Section we will prove in the latter norms the convergence of the sequences C n , F n . The limiting 
covariances will thus turn out to be in C p , V» > 0. 

We recall the definition of the lattice derivative V £ in (|2.7j) . In particular S = {ex, ■■■e-d\ is the 
standard basis of unit vectors. For e£5, Ve is the forward partial derivative and for — e £ S, Ve 
is the backward partial derivative. We define the n th lattice derivative 



V£ e =V ei V e 



Let X be a connected open set in H d . We define 

x £ = xn {eZ) d 

We now define the lattice Sobolev norm ||.||H fc (x e ) of a function / by 



k 

2 _ 

i=0 ±^,...,±^€5 



Lemma 5.1. Zei g be a C°°(R d ) function. Then for every k > i/iene exists a constant Ck 
independent of e such that 

\\fg\\H h (X e ) < Cfclbllc fc (R d )ll/llH fc (x £ ) 
For every k > 0, and any Li function g, 

ll/*2lltf fc (( £ Z)<i)) < H^llL 1 (( e z)««)ll/llH fc (( e z)<')) 
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Proof. To prove this we take the square of the norm on the left hand side and then use the (eZ) d 
lattice modification of the Leibniz rule : 

V £) e(/<?) = (V £ J)g + fV £ ,e9 + eV £ JV £ ,eg (5.2) 

Derivatives on g are bounded in the C k (R, d ) norm and all e dependent constants can be ma- 
jorised by setting e = 1. This proves the first inequality. The second inequality, which is a form of 
Young's convolution inequality, is proved exactly as in the continuum. □ 

In the following we will exploit a lattice version of elliptic regularity. Let / be a bounded 
function in (eZ) rf . Let U(R) C R d be an open cube centered at the origin and of edge length R. 
Let U £ (R) = U(R) n (eZ) d . Let O C U(\R) be an open connected set. Then Q £ C U £ (\R). Define 

h a (x)=V^(x,f) 

Recall from Section 2 (see (|3.3JI . (|H.4|) ) et seq.) that h a (x) is the solution of the Dirichlet problem 

(-A £ + a)h a (x) =0: x € U e 

h a {u) = f(u) : u G dU £ (5.3) 
and that the maximum principle holds, because a > 0. We have 

Proposition 5.2 (lattice elliptic regularity). With Q £ defined as above, andV k > 

\\h a \\ Hk (n £ ) < C R (l + a)-^\\f\\ L ^ ii£Z)d) . 

Remark : This is well known in the continuum. For completeness we give a proof of the lattice 
version in Appendix [SI 

We choose e = e n and apply the proposition to 



(^ njm=1 (L°)/)(x') := / dzc^Mx'-^n: (lA x 'if) 
which is the j = n term in the product (|3.28|) defining A^- 

Corollary 5.3. Let Q £n be as in Proposition^^ with R = 1. Then for every k > and every 
n>\ there exists a constant C^^l independent of e n such that 

\\A a entm=1 (L°)f\\ Hk{QEn) < C ktL (l + a)-h\\f\\ Leo({enZ)i) . (5.4) 

l|V| 1 ,...,e fc ^, m= i(^ )/llL oo (( £ „zn < C kjL {l + ayHf\\ L ^ {enm . (5.5) 

Proof. The integral over z in (j5.4l3.9f) can be restricted to U £n (2) because the range of g\ is 1/4. 
Therefore 

\\ A £ n , m =i( L °)(-J)\\H k (n £n ) < 

/ dzc en-i\\9l\\cHm\\K: n (l,z)^^^ H ^ n ) 
JU en {2) 

< c ktL (l + a)-^\\f\\ LX{{enZ)d) 
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We have used the first inequality of Lemma 15. 11 absorbing the C k (R d ) norm of g\ in the constant 
since g has been fixed once for all , and then using Proposition 15.21 

To prove (|5.5|) : By the embedding of high degree Sobolev space into L^, (Lemma lB.l() . reviewed 
in Appendix El we pass from (|5,4j) to 



noting that (|5.4j) applies to any translate x + Q £n of Q £n . 



(5.6) 
□ 



Lemma 5.4. For every integer k > 1, and every n > 1, ] a constant ct } L independent of e n such 
that , 

\A^ m=1 (L )(p)\<c k , L (l + a)-h(-A £n (p) + l)- k (5.7) 

Proof. This is the essentially the standard proof that the Fourier transform of a smooth function 
of compact support has rapid decay. Let 



h(x,p) 



Al im=1 (L°)(x,du)e-^ 



Then 



-ipu 



(-A e » + l) K h(x,p) = / A a en:m=1 (L°)(x,du){-A Ut£n + 1) V 

J(e n Z) n 

where the u subscript on A Uj£n indicates the variable it differentiates. Since A Uj£n is a linear combi- 
nation of lattice translations under which the lattice is invariant, the exact analogue of integration 
by parts is valid and we continue with 



/ {-A u , £n + l) k A^ m=1 (L°)(x,du)e- 

J(e„Z) n 



i pa 



By the translation invariance (see (jH.K)j) ) of A® n m=l (L°)(x, du) we can change the derivatives to 



/ (-A XjEn + l) k A a £n!m=1 (L°)(x,du)e 



-ipu 



'(e„Z) 

By ()5.5() in Corollary 15.31 with f{u) = exp(ip.u), 



'(e n Z) 

Collecting these relations we have 



/ (-A se>en + l) fc ^ n)rn=1 (L°)( a ;,d«)e- 

He n ZY 



ipu 



<c k:L (l + a) 2 



<c fcjL (l + a) 2 



-A £n (p) + l) K h(x,p) 

By setting x = we finish the proof. 

Fourier transforms are naturally defined on the Brillouin zone 

B £ = [- 7r /e,Tr/e] d 
There is a constant c independent of e such that 

p 2 > — A £ (p) > cp 2 for p E B £ 
which follows from t 2 /2 > 1 — cost > ct 2 on [— it, ir]. 



□ 

(5.8) 
(5.9) 
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Theorem 5.5. Vn > and \/k > 0, 3 a constant Ck,L independent of n such that 

K(p)\<c k , L (l + a)-\l+p 2 )- 2k for P eB en (5.10) 

\\K\\H k ((e n zr) <c KL (l + a)- 1 (5.11) 

Proof. iFrom (l2~TU)) and (ETUI) . 



c 



0<^»<(-A £ J-<^ 

Combining this with ()3.16)) and the continuity assertion of Lemma 13.11 we have 

f a En (p)<c L (l+p 2 r l 

Case n = 0: It is sufficient to prove that T £Q (p) is bounded by C(l + «) _1 with C uniform in p 
because B eo is bounded. Referring to 1)3.1612. 1U|) we find that 

if* h)\ < 2 |1 - i °<> )1 < J 1 -KM j4 a (Q)-4 a o(p)l 

£oWI " a-A £0 (p) " « a-A £0 (p) 

The first term is continuous at a = by Corollary 14.21 Therefore it is bounded by C(l + a) . 
The second term is bounded by C(l + a) -1 using the same argument (existence of moments of the 
Poissson measure) as in the proof of Lemma 13. 11 

Case n > 1: by Lemma 15. 4) and the bound |A £n m (p)\ < 1 which we use for m > 2 

K(p)i 2 <c fc , L (i+arv + ir fc 

Using these estimates in ()3.29|) we obtain 1)5.1 Oil . 

Proof of (|5.11|) . We have the easily established bound 

\P 2 + A £n (p)\ <o(i)4H 4 

so that 

< 1 - A £ » < (1 +P 2 )(1 + 0(1)4/) 

and hence for any m > 

(i - A £ »nf»i 2 < o(i)(i +/) 2m if»(p)i 2 < Cfc , 3 , L (i + oj-^i +p 2 )- ? 

for any g > by choosing k in ()5.10)) sufficiently large. Taking q > d proves the theorem because 

□ 

Now turn to the Levy fluctuation covariance given in 1)4.3)) . Using the bounds provided in 
Theorem 15.51 we get 

Corollary 5.6. For < a < 2, all k = 0,1, ... , and aZZ n > 0, 

ll r n|lH fc ((e n Z) d ) < c fc,L (5-12) 
where the constant on the right hand side is independent of n. 
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6 Convergence 

Theorem 15 . 51 and Corollary 15 . 61 provide uniform bounds in Sobolev norms for fluctuation and block 
covariances. In particular they are uniform in the lattice spacing e n . We will now prove that these 
sequences converge to their formal continuum limits. Continuum objects have the subscript c in 
place of e. Thus, as in (|2.9|) . 



G a c {x-y)= I * e*<*-*)(9»(p) (6.1 



Recall from Section |3] that 



d d p 
Rd W 

G a c (p) = (a + p 2 )- 1 (6-2) 



u(R)^u c (R) = (-^,^) d cn d 



represents an open cube of edge length R. In analogy to (|3.23|) with c c := 1 we define the continuum 

average A a {R m )(x, , du) by 



L 



duA a {Rm){x,du) f{u)= I dzc c g m {x- z)V a v(R -\(x,f) (6.3) 



where 



n c (R, z) &f)= [ duV a Uc{ ^ z) {x } u)f{u) 
JdU c (R,z) 



is the solution h"(x) to the continuum Dirichlet problem 

(-A c + a)h%(x) = : x G U C {R) 



h a c {x) = f{x) : x £ dU c (R) (6.4) 
With these notations, the Fourier transform of the continuum analogue of (|3.13l) is 

f«(p) = G a c ( P ) - \Al (R )(p)\ 2 G a c ( P ) (6.5) 

and that of (EHHl) is 

n 

fc,n= II \A a c , m (Rm)\ 2 r a c (6.6) 
m=l 

The Lemmas, Propositions, Theorems and their Corollaries of sections 3,4 and 5 remain true in the 
continuum with the following caveat : in the continuum the uniform Sobolev bounds of Theorem 
15.51 and Corollary 15.61 hold only for n > 1 . Aside from this caveat their proofs are identical and 
need no repetition. When referring to them for the continuum objects we shall simply mention 
them as the continuum analogues of the relevant results for the lattice. 

Recall that B £n is the Brillouin zone, defined in l|5.8j) . The main result is 

Theorem 6.1. For every integer k > 0, 

K, n - r£. (6.7) 
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in Hk(R, d ). Moreover, for every fixed lattice (eiZ) d , < I < n the restriction ofY^(x) to {e{Zi) d 
converges to the continuum Y^ st (x) restricted to (eiZ) in the Sobolev norm 

\K* ~ r nlk(( £i zn - as n - oo (6.8) 

Moreover multiple lattice derivatives of Y® converge to the corresponding continuum derivatives of 
in the Loo((£;Z) d ) norm. 

For the Levy finite range decomposition (|4.2j) we apply the last theorem to (|4.3[) and obtain 

Corollary 6.2. For all k>0, 

||r C) * - r n ||^ fe (( £;Z)d) -> as n -> oo (6.9) 

Moreover multiple lattice derivatives of Y n converge to the corresponding continuum derivatives of 
Y Cj * in the Loo((e;Z) d ) norm. 

We now give some Lemmas which will be employed in the proof of Theorem 16.11 In the 
following lemmas we consider continuum functions / : H d — > R and use the same symbol for the 
lattice function / : Z d — ► R defined by restriction. Continuum and lattice integration are to be 
distinguished by the domain of integration. 

Lemma 6.3. 

\\f(-+h)-f(.)\\ Hk{nd) <c\h\.\\n Hk+l{ n^ 

where \h\ is the norm in Zi d C YL d . Moreover if \h\ < R we have 

||/(. + h)- f{.)\\ Hk{ U c{ R)) < ^1^1.11/11^(^(2^)) (6-10) 

Proof. See Theorem 3.3 on page 42 of |Agm65| . 

□ 

Lemma 6.4. Define 

Q = -A e - (-Ac) 

Then we have for every k > 0, 

HQ/lk(R«) < ce \\f\\ Hk+3m (6.H) 
where the constant c is independent of e. 

Proof. Ve, e is the forward lattice derivative in (eZ) d , and V* £ , its L 2 ((eZ) d ) adjoint, the backward 
lattice derivative. Forward and backward derivatives commute. V e . c is the continuum derivative in 
direction e and the adjoint is V* c = — V e c . A calculation shows that 



Therefore 



V* £ V e>£ /(x) = / dt f ds (-Vl jC ) f(x + [t- s]se) (6.12) 
Jo Jo 

Qf(x) = - J2j dt I ds ( v lcf( x + [* - S M - v|, c /(x)j 
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since 

-A £ = ^V* i£ Ve, e 

e 

the lemma is proved by taking norms under the integrals and using Lemma 16.31 (which gives the 
factor e). 

□ 

We can now describe the main idea. Let / : dU c (R) — > R be a continuum function / : R d — * R 
restricted to the continuum boundary of U C (R). We need to estimate the difference between the 
solution h". to the continuum Dirichlet problem and the solution h® to the lattice Dirichlet problem. 
This will be done by restricting h® to the lattice. The restriction also solves a lattice Dirichlet 
problem, but with a non-zero right hand side involving Qh^, which by Lemma 16.41 is 0(e n ). Thus 
we consider 

h a £n {x) - K(x) : x e iuf ) C U e (*) (6.13) 
h a en (x)-h a c (x)=0 : xGdU En (R) 

in Sobolev norms. The difference satisfies zero boundary conditions because we will arrange that 
the lattice boundary points dU £n {R) all lie on the continuum boundary dU c (R) and both solutions 
have boundary values / restricted to dU £n (R). 
The lattice cube 

U £ (R) = U(R) n (eZ) d 

has as its boundary 

8U e (R) = {y g U £ {R) : \x-y\=e, some x G U £ (R)} 

iFrom Section |21 we have e n = L~ n , L = 2 P , p > 1. In addition we now choose , in accord with 
= R m = L^™- 1 ) : < m < n. Then not only do we have 

U £n (R)cU c (R) :=U(R) 

but also 

dU £n (R) C dU c (R) (6.14) 

This last statement follows from the observation that e n = 2~ np and 4p = 2 _mp+p_1 so that for 
< m < n we have ^ £ e n Z C R. This means that for d — 1 the boundary points of U £n (R m ^) 
coincide with the boundary points of U{R m ). For d > 1 we are in cubes and the above reasoning 
easily generalises to ()6.14[) . 

Lemma 6.5. Let 

f : K d —> R, 

Then for 

R = L-t™- 1 ) : < m < n 

we have 

ll V e 1 ,...,e fe ( /i e„ ~ h c)\\ LaD (U en (f)) ^ c L,R,k £ n II / II L°° (dU c (R)) 

where the constant ci t R t k is independent of n. 
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Proof. See Appendix A. 

To proceed further we need a formula. First define a new finite difference derivative that acts 
on continuum functions by : 



Ve,ef(z) = [ eft (1 - t)V e . c f(z + te) 
Jo 



Then we have 
Lemma 6.6. 

d-l 



[ dz f[z) - [ dz f(z) = £ V / dz V e c f(z) (6.15) 

Jn. d J{eZ) d ~q J(eZ) d -J'xRJ 

Proof. We obtain case d = 1 by 

/ /(*) - f dz f(z) = V / dx /(*) - / /(z) 

dz a' 1 f dx {f{x) - f(z)) =e f dz V e , £ f(z) 

eZ J[z,z+e] J eZ 

and then the general case is obtained by applying this formula iteratively. □ 
The Fourier transform of A% m (R m )(0, du) is given by 



where 



At n , m (Rm)(p) =c £n _ m / dzg m (z)e w ' z h a {z,p) 

J(e n Z) d 



K n , m (z,p)= [ V a UEn(Rm) (z,du)e^ 
JdU £n (R m 



Likewise, there is the continuum Fourier transform A^ m (R m )(p) defined by the same formulas with 
e n ,e n - m replaced by c and with c c := 1. 

We wish to estimate the difference A^ m+1 (R m+ i)(p) — A^ n ^ m+1 (R m+ i)(p). This is provided by 
the following Lemma : 

Lemma 6.7. For all integers k > 0, and n > 1, 3 a constant Ck,L,m independent of n such that 

\Ac, m+1 (R m+ i)(p) - A a E ^ m+1 (R m+1 )(p)\ < c k)Ltm e n 
Proof. It is easy to see using Lemma 16.61 that we can write 

A a Ctm+1 (R m+1 )(p) - A a £nim+1 (R m+1 )(p) = R 1 (p) + R 2 {p) + R 3 (p) (6.16) 

where 



Rl(p) = (1 - c £n _ 1 _ m ) / dz e ip - z g m+1 {z)h a +l {z,p) 

J(e„Z) d 

#2(p) = c e«-i-m / dz e- ip - z g rn+l {z)(h a +l {z,p) - h a en +1 {z,p)) 

J(e n Z) d 
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d-1 . 

Rzip) = e n V / dz e~ ip - z Ve c (g m+1 (z)hl m+1 (z,p)) 

~7 J(enZ) d -jxW 

We observe that 

|1 — C r . I < CTjEn 

as follows from the definition (|3.7|) . Lemma 16.61 and the fact that g is a smooth function in R rf of 
compact support. The integral is bounded by 0(1). 

For the term involving R2 we use Lemma 16.51 with f(u) = exp(ip.n) which produces the small 
factor e n in the bound. 

Finally in R3 the 0(e n ) factor is already there. In the integrand h^ m+1 (z,p) and derivatives are 
Lqo bounded on the support of g m +i by (the continuum versions of) Proposition ^, 2l and Lemma lB.ll 
with /(it) = exp(ip.it). 

□ 

Proof. (Theorem 16. 1[) Let e > and fix any p ^ 0. Fix any a > 0. We will first prove that 

|f^(p)-f»(p)|<e 

for all sufficiently large n. 

Recall that A a en m (p) and A^ nm (p) are Fourier transforms of defective probability measures 
supported in a cube of side R m . Now 

\l-A a £n>m (p)\ <|l-/ A a enim (0,du)\ + f ^ nim (0,d«)|l-exp(-ip.x)|. 

Therefore, by |1 — exp(— ip.x)\ < c|p.x|, and Corollary 14.21 

|1 - A a Enjm (p)\ < cR m \p\ + caR 2 m (6.17) 

Note that the same bound holds in the continuum, because Corollary 14.21 remains true in the 
continuum. 

^From the definition of n we get 

r c ,n+l — ^c,n = (j A:,n+l(^n+l)| — ^j^c,n 

whence, using the continuum analogue of (|6,17l) . 

Kn+1 " fc,nl < (2cRn+l\p\ + 2 C ai^ +1 ) |f £ n | 

Now r° n satifies the uniform bound of Theorem 15. 5[ R n +i decreases geometrically with increasing 
n. Therefore we see from the previous inequality that T^ n form a Cauchy sequence. This proves 
the existence of the limit T® ^ satisfying the bound of Theorem 15 . 51 and the first part of Theorem 16. II 
has been proved. 

Choose n sufficiently large so that 

\%&)-%m<\ ( 6 - 18 ) 
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f 1 " (p) is uniformly continuous in p by Lemma 13.11 Therefore, as shown in the course of proving 
Theorem 15.51 . there exists a constant 7 independent of e n such that 

|f£»l<7 (6-19) 

Now (|6.17j) implies that 

|l^£„,ml ~~ l^c,ml| — 2ci? m |p| + 2caR m 

Since R m — > geometrically fast and |A" m (p)| < 1, we can choose N depending on e such that for 
all n > N, 

n n 

1 n n Kmi a i<i (6-20) 

m>AT m>A r ' 

By Lemma 16 . 71 there is a constant Cat such that 

N N 

I II \ A lm\ 2 - II l^-mH < ^ < i ( 6 - 21 ) 
m=l m=l 

for all sufficiently large n. Finally, 

\t a en -T a c \<- A (6.22) 

by the definitions (|3.13|) , the explicit Fourier transforms 12 . 1 Ull6 . 2|) and Lemma 16.71 

^From the definition (|3.29|) we see that the four inequalities 1)6.191 16.201 16.211 16.22)1 imply that 

Kn-K\<J (6-23) 
()6.23j) and ()6.18)) establish the pointwise convergence 

Kip) - K,M (6-24) 

By the dominated convergence theorem using Theorem 15.51 for domination, we have, for any 
fixed compact set X C R d in momentum space and any k 



[ dp\flJp)-f a n (p)\ 2 (l+p 2 ) k ^0 
Jx 



as n — > 00. This proves 1)6.8)1 of the theorem, because we can choose X to be a fixed Brillouin zone 
in the dual of (e;Z) d . The convergence in the -Loo((£zZ) d ) norm follows by Sobolev embedding (see 
Appendix B). 

□ 
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A Lattice Elliptic Regularity 

Suppose that h solves 

(a-A £ )h = g (A.l) 
on (eZ) d and ip has compact support. Then, on (sZ) d , 

(ph(a — A e )h = iphg 

Integrate over (eZ) rf . By the definition of — A £ this can be rewritten as 

a / dz (ph 2 + / dz y^(V e <ph)(V e h) = / dz iphg (A.2) 

Surprisingly, calculation shows that this can be rewritten as 

a dz <ph 2 + \ V / dz <p(X7 e h) 2 = / dz [phg + \(A e ip)h 2 ] (A.3) 

JW 2 ^ s J(ezy J(ezy 2 

In other words with forward and backward derivatives the lattice gives the same formula as the 
continuum without corrections that go to zero with e. For a > and for ip > 0, (|A.3|) implies 

l J2 f dz p(V e h) 2 < f dz [phg + l(A £ p)h 2 ] (A.4) 
I dziph 2 <-[ dz [phg + -(A £ ip)h 2 ] (A.5) 

J{eZ) d a J{eZ) d 2 

Let 

\K,i = 2 ~ j E / dz\V{ h\ 2 p 

By applying j finite difference derivatives to (jA.lj) we find that (|A.4|) is also true for h and g 
replaced by derivatives of h and g. By the Cauchy- Schwartz inequality on the first term in (|A.4|) . 

K d +i<[h] v M™ + [ h ]\A eV j 

Simplify the first term using the inequality ab < |(a 2 + b 2 ) and use the resulting inequality to 
iteratively reduce the order of the top derivative in the Sobolev norm, 

\\h\\l, h --=J2 2 ~ j E / d *\V j ei ,..,eM 2 <P (A.6) 

We obtain, 

INlL^IHlLo + NI^-i ( A - 7 ) 

where <po is the final member of a sequence ipj of non-negative functions chosen such that 

1 1 . , . 

(Pj-i > gVj + 2 Ae ^'' Vj-l > fj ( A -8) 
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Proof. ( Proposition 15.2)) . We are given that h solves (jA.lj) in U £ (R) with g = 0. We can estimate 
the L2 norm in U £ (R/2) in two different ways. Firstly, For all a > we can use the maximum 
principle \h\ < \\f\\ L oo^ gu ^ R)) 

\\h\\L 2 (Ue(R/2)) < cRd/2 \\f\\L°°(dU £ (R)) 

Secondly, for a > 1, we can choose ip in (TOl) to be one on U(R/2) and zero outside U(R) to obtain 

IHlLoW*)) < CR*~ l \\h\\l, {Uem ^ A - 9 ) 
and then use the maximum principle to bound the right hand side by the L°°(dU £ (R)) norm. 
Therefore, for any smooth po supported in U e (R/2) 

\\h\\l 0! o<c R (i + ay^f\\l^ dUem (A.10) 

Let (fk, <Pk-i, ■ ■ ■ j Po ^ be C°° continuum functions with compact support in U(R/2) such that 
<Pk = 1 on and l)A.8j) holds. Apply (|A.7j) and ()A.10j) to obtain Proposition E21 □ 

Remark on Exponential Decay: The correct a dependence is exp(—0(^/aR)). We outline how 
to do this using a method suggested by |Agm82| . (|A.5|) . for g = 0, can be rewritten as 

dz wh 2 < 

(ezy 

w:=(a- 2 A eV 

Consider the choice 99 = exp(— u) with 



u(x) = —\fa~v 1 + x 2 
V ) 2 v V 

Let Ao be the continuum Laplacian. The finite difference Laplacian A £ (p can be written as an 
integral over Ao^: for example, in one dimension 

-1 r l 



A £ ip(x) = ds dt Aqp(x + [t — s]e) 

Jo Jo 



Using this we find A £ u > and w > c(a,e)ip with c(a,e) » exp(— y/ae). Now replace exp(— u) by 
<p = exp(— u)ip where ip is a smooth, positive, monotonic decreasing function such that ip = 1 on 
and vanishes outside U(R). Then w > c{a,e)p > on C/(|i?) so by taking the part of the 
integral where w to the right hand side of the bound and discarding part of the integral where 
w > we get 

c(a,e) dz ph 2 < dz {-w)h 2 < \\ff L ^ idUsiR)) / dz \w\ 

Ju e (R/2) Ju e (R)\U E (iR) Ju £ (R)\U e (iR:) 

which gives decay because 

<f\u(R/2) > ex P(0(VaR))p\ u{R) \ u{ 2 R) 
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Preparation for proof of Lemma \6.5\ Suppose (|A.1|) holds in a domain U £ (R) and h vanishes 
on dU £ (R). Then, with 93 = 1, a > 0, l)A.2j) becomes, 



V / (Veh)(V e h) < dzhg (A.ll) 



>U e (R) JU e {R) 

We estimate the right hand side by the Cauchy- Schwartz inequality and substitute the result into 
the Poincare inequality, which is 



\L*(U C (R)) 



<C R Y] [ dz (VJi)(VJi) (A.12) 
Ju, 



where Cr is independant of e. Then 

\\ h \\L 2 (u e (R)) < Cilllf IU 2 (C/ e (fl)) (A.13) 
Returning to l)A.7|) and using (jA. 13|) we obtain 

Lemma A.l. Let <p > be a C°° continuum function with compact support in U(R). By restriction 
it defines functions on all lattices, denoted by the same letter. Then for e sufficiently small, there 
exists a constant Cr )V such that the solution to \A.1\) . with zero boundary conditions on dU £ {R), 
satisfies 

\\h\\ip,k < CR^WgWH^xiUeiR)) 
where the constant Cr^ is independant of e. 

Proof. (Lemma 16.5(1 Write as in lemma 4-3, —A £n+1 = — A £n + Q. Let h = h® n — h* 1 . This has 
zero boundary conditions on dU £n (R) and 

(-A £n + a)h(x) = g{x) : x G U £n (R) 

where g = Qh^ n+1 . Now apply Lemma IA.1I and estimate g by Lemma 16.41 followed by Propo- 
sition 15.21 The proof is completed by Sobolev embedding, see Lemma IB. II below, taking k in 
Lemma lA . II sufficiently large. □ 

B Sobolev Spaces on the Lattice 

Lemma B.l. Let L d := [0, l} d , L d = L d n (eZ) d and X be any subset of L d . Then 

\\u\\l°°(x) < C d \\u\\ Hd ^d-) 
Proof. A smooth continuum function u(x), where x = (x±, . . . , Xd), satisfies 

xi... x d u(x) = dyi... dy d di... d d (yi . . . y d u{y 1 , y d )) 
Jo Jo 

and therefore, for x £ L d := [0, l] d , 

" \x 1 Cd .x d \ MHd{Id ~> 
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The same proof adapts to the lattice (eZ) d with integrals and derivatives being replaced by sums 
and finite differences so that for u : if — > R, 

sup \u(x)\ < C d \\u\\ Hd{Ii) 

x£j £ 

where J £ := [1/2, 1]^. For d = 2, I 2 is the union of [1/2, l] 2 , [0,l/2] 2 , [1/2,1] x [0,1/2] and 
[0, 1/2] x [1/2, 1] and by symmetry the same bound holds with J £ replaced by any of these boxes. 
The same argument applies for all dimensions d, which implies the lemma. □ 

Lemma B.2 (Poincare Inequality). Let u : if — ► R be any function vanishing on the boundary 
of if. There exists C independent of e such that 



L 



dx\u\ 2 < CHViill^dj 

Proof. Let Ui(x) be the finite difference partial derivative with respect to component Xj of x. 
Then, for i = 1, . . . , d, 

\u(xi, . . . ,x d )\ < \ dxi\ui(xi, . . . ,x d )\ 

lie 

The right hand side is a function of all components of x except Xj. Take the product over i followed 
by 2/d root or power. 



2 1 

d 



Hx)\ 2 <n(/ ^ki) d <n(/ dxi\ Ui \ 2 ^ 

Integrate both sides over (xi, . . . , x d ) € if and use the Holder inequality on the right hand side. 
/ cfe|-u| 2 < TT ( / dx I dxi\ui\ 2 ) < / dx"S^ \ui\ 2 

J If i\JIf Jle J Jli i 

because the extra integral integrates to unity. □ 
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